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Abstract 

An attempt is made to clarify the ballistic non-linear sigma model formalism recently pro- 
posed for quantum chaotic systems, by looking at the spectral determinant Z(s) = Det(l — sU) 
for quantized maps U G U(iV), and studying the correlator u>u(s) = J dO \Z(e ie s)\ 2 . By identi- 
fying U(iV) as one member of a dual pair acting in the spinor representation of Spin(4iV), the 
expansion of u>u(s) in powers of s 2 is shown to be a decomposition into irreducible characters 
of ~U(N). In close analogy with the ballistic non-linear sigma model, a coherent -state integral 
representation of uu(s) is developed. For generic U this integral has ( 2 ^) saddle points and the 
leading-order saddle-point approximation turns out to reproduce iou(s) exactly, up to a constant 
factor. This miracle is explained by interpreting u>u(s) as a character of U(2iV), and arguing 
that the leading-order saddle-point result corresponds to the Weyl character formula. Unfor- 
tunately, the Weyl decomposition behaves non-smoothly in the semiclassical limit N — > oo, and 
to make further progress some additional averaging needs to be introduced. Several schemes are 
investigated, including averaging over basis states and an "isotropic" average. The saddle-point 
approximation applied in conjunction with these schemes is demonstrated to give incorrect re- 
sults in general, one notable exception being a semiclassical averaging scheme, for which all loop 
corrections vanish identically. As a side product of the dual pair decomposition with isotropic 
averaging, the crossover between the Poisson and CUE limits is obtained. 



1 Introduction 

One of the striking characteristics of a quantized chaotic Hamiltonian system is found in 
the correlations inherent in its spectrum at small energy differences. Extensive numerical 
work has shown that various quantities (such as the nearest-neighbour spacing distribu- 
tion and the two-level correlation function) of a quantum chaotic system are universal: 
their behaviour coincides with that of a Wigner-Dyson random matrix ensemble of the 
appropriate symmetry class 0. This property, first noticed in billiards, was found to ap- 
ply to many chaotic systems, including symplectic maps. In contrast, if the dynamics is 
integrable (in the sense that the 2 /-dimensional phase space foliates into /-dimensional 
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submanifolds invariant under the Hamiltonian flow) the generic behaviour of the eigenval- 
ues is expected to be that of independent random variables, so that their correlations 
are in the Poisson universality class. 

The present paper will be concerned with quantum maps, i.e. with quantizations of 
some canonical transformation cf) : M — > M of a compact symplectic manifold M. We 
assume that the problem of quantization itself has been tackled, so the phase space has 
been prequantized into a Hilbert space Ti^ of dimension N ~ fir 1 , and the quantum map 
acts on it as a unitary operator [[33|, |3|. With respect to a basis of TCn this operator is 



represented by an N x N unitary matrix U^. The latter has a semi classical limit, in the 
sense that traces of its powers can be estimated in terms of classical periodic points 
For a system with one degree of freedom, the Gutzwiller-Tabor trace formula reads 

Tr(U^ N ) 7V ~ 00 N dim ^/ 2 A p e iN ^ , (1) 

pCFix(</>") 

where p is one component of the set of n-periodic points; for an Anosov system, it is an 
isolated point (dim(p) = 0), whereas if the dynamics conserves energy, p is 1-dimensional. 
<3> p and A p are purely classical quantities related to the dynamics around the set p. 

The quantum spectrum consists of the A" eigenvalues (pseudo-energies) {e 1 ^},^^...^ 
of UfaN- The first analytical estimates of the two-level correlation function (which is the 
Fourier transform of the form factor F{n) = \Tr(U2 N )\ 2 ) for such spectra were based 
on the above trace formula, combined with some known ergodic properties of long peri- 
odic orbits 0. In the present paper we focus attention on another statistic, namely the 
autocorrelation function of the spectral determinant: 



fV( 7 ) = T N/2 !^ — Det(l - 7 e^[/)Det(l - e~^lP) . (2) 
Jo 2vr 

(The parameter 7 will be a complex number close to unity, with the scaling I7— 1| ~ 
This correlation function has already been considered || [7]] for chaotic versus integrable 
quantum maps, and the same universality was observed as for the form factor or the 
nearest-neighbour distribution. A semiclassical analysis of this correlation function was 
performed using the Gutzwiller trace formula in f/J. 

The computation of correlation functions from the trace formula ([I]) always requires 
some sort of averaging. In the semiclassical theory of the form factor F(n) one wants to 
use the so-called diagonal approximation, neglecting the off-diagonal terms in the double- 
sum over periodic orbits [f|. To justify this step one must average over energy or some 
family of systems: one needs slight variations in the classical actions $ p to make the 
phase interferences ~ e lAr< -* p_<I> p'- ) average to zero. The need for averaging was emphasized 
in J6|, [8| , where it was pointed out that the spectrum of an individual quantum system is 
too noisy to allow universality to be seen in its bare form factor. In addition to the noise 
problem, there exist some quantum chaotic systems with arithmetic symmetries, which 
lead to periodic orbit degeneracies and non-universal spectral correlations || [10| . (Such 
systems are non-generic, however, in any decent space of smooth maps.) 



2 



Thus universal behaviour is expected only in the generic case, and to make a correct 
mathematical statement about universality of the spectral correlations of a general system 
one ought to define the precise meaning of the word "generic" . (In the case of integrable 
systems, the spectral correlations could sometimes be studied directly, by utilizing the 
explicit expressions for the eigenvalues; two-point correlations were shown to be Poisson 
for a rather subtle set of parameters Q.) We can avoid the issue of genericity by averaging 
the correlation function over some set of quantum maps. That is, we specify a measure 
(1Vn(U) on the unitary group U(iV), and the function to be studied then reads 

(^(7)) = / dV N (U) Qu( 7 ) • 

JV(N) 

We want this measure to be very concentrated (or "local" ) around the quantum map U^n 
(see Sections [4.3| , §). In the course of this article, we will also consider cases for which this 
measure has a broader support (Section ||), including the case where the measure is the 
Haar measure on U(iV). We are then dealing with the circular unitary ensemble (CUE), 
for which the determinant correlation function has been thoroughly investigated |6], |Tl| . 

The Gutzwiller trace formula has the attractive feature of relating quantum to classical 
properties, but its use for estimating the spectral correlation functions still raises ques- 
tions. The problem is that the formula is rigorous in general only for times shorter than 
the Ehrenfest time: n < logiV (and iV — ► oo). Yet, long-time traces (n ~ N) are needed 
to obtain spectral correlation functions at the scale of the mean level spacing where uni- 
versality emerges. The diagonal approximation, which assumes statistical independence 
of the different periodic orbits, is unsatisfactory at large times where the exponential 
proliferation of periodic orbits clashes with the finite (N) number of eigenvalues: the 
classical information is then overcomplete, which implies some sort of hidden correlation 
between the contributions from classical orbits. A recipe to overcome this difficulty has 
been devised by Bogomolnyi and Keating Hl5| , but so far lacks rigorous justification. 

To bypass these problems, a second approach to estimate spectral correlations has 
recently emerged, inspired by the study of disordered metals. It consists in expressing 
the correlation function as a quantum field theory (or functional integral) of the type of 
a non-linear sigma model (NLcrM). One then tries to analyze the functional integral by 
standard field-theoretic methods such as perturbation expansion, saddle-point analysis 
and the renormalization group. This approach was first applied successfully to systems 
with disorder, where the dynamics is governed by a diffusion operator Hl2| . The formalism 



was later extended to the "ballistic case" [|13|, [L4[] , and quantum correlation functions were 
put in relation with the spectrum of the Frobenius-Perron operator (i.e. the evolution 
operator for classical densities). Although quite elegant, this approach suffered from 
several drawbacks. Among these are the appearance of unwanted zero modes around the 
main saddle point, and the problem of "mode locking" [p^] . Besides, the results do not 
exactly agree with the correlations calculated numerically for the Riemann zeta function 
(the prototype of a quantum chaotic spectral determinant) |15| ; nor do they explain the 
non-generic spectral correlations featured by systems with arithmetic symmetries. More 
recent treatments of the ballistic NLaM have also stressed the need for averaging over a 
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smooth disorder if one wants to avoid the above problems [fTo] , |T7| . 

In an attempt to resolve these uncertainties, we have adapted the latter approach, 
which had originally been conceived for Hamiltonian systems, to the case of quantum 
maps U<j, t N- Our objective was to prove the universality of the determinant correlation 
function fit/ (7) (Section |2.2| ) upon averaging w.r.t. a suitable measure on U(iV). This 
correlation function is easier to treat than the two-level correlation function, as it does 



not require the use of a supersymmetric representation |14| but can be expressed as 
an ordinary c-number integral over a finite-dimensional manifold Aijy (Section |2.4|) . We 
write this integral in the form 

nu(l)= I dQe-sWV, (3) 

J Mn 

where 5(7, U, •) is called the effective action. 

For the reasons stated, we will consider averages of Qu with respect to certain prob- 
ability measures AVn(U). The averaged correlation function, denoted by (^^(7)), can 
still be obtained by integrating the Boltzmann weight given by an effective action: 

(Q u ( rf ))= I dQe- s ™^' Q) . 

J Mn 



To estimate these integrals, we apply the same technique that was used in [pj: we expand 



Sav(7, Q) up to quadratic order around its saddle points Q C rit, and perform the Gaussian 
integrals. The result obtained in this way, 

<^))|,p. exp. = E ( Det ^(7, Qcrit)y 1/2 e-^«-) , (4) 

Qcrit 

is called the leading-order saddle-point expansion of the integral. 

Owing to the absence of a large parameter in front of the action S aY , the expansion is a 
priori not justified mathematically. A more careful treatment should in principle include 
perturbative corrections around each critical point (we actually compute the expansion 



up to two loops in a particular case, see Section [4~4] ). 

We have succeeded in computing the leading-order term for a few averaging schemes. 
For an individual matrix Un we can actually reproduce the exact value of the correlation 
function (|j) in this way (Section §). In Section [4.3| we define a "semiclassical" averaging 



scheme, which we think is a good candidate to obtain universality of correlations [ 21 
unfortunately, in that case we can only compute the contributions from the two standard 
saddle points. 

In order to test the leading-order saddle-point approximation, we selected a sequence 
of statistical ensembles (i.e. a sequence of measures &Vn{U)) for which the averaged 
correlation function can be computed exactly, and compared the exact result with the 
saddle-point approximation for the corresponding effective action. All these ensembles 
are U(iV)-rotation invariant, that is, we first average over all bases of 1~Ln (Section |5|), 



then possibly over the spectrum of Un (Sections In most cases, the saddle-point 
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expansion of these ensembles yields erroneous results. We still hope that the expansion 
is better behaved in the case of local averages, like the semiclassical one. 

These disappointing results seem to challenge the use of NLcrM methods for the study 
of quantum ballistic systems, unless our understanding and control of these methods 
significantly improves. In Section || we introduce a U(iV)-isotropic local averaging scheme 
which we treat by an alternative method; unfortunately, this scheme does not discriminate 
between the different universal behaviours that are expected for chaotic versus integrable 
maps. Nevertheless, we use it in Section ^T] to compute the correlations along a crossover 
between the Poisson and CUE universality classes. 



2 Algebraic manipulation of Vtjj 
2.1 Fourier decomposition of Qjj 

We first remind the reader of some known results concerning the correlation function Qu 
|| [7]. The spectral determinant of U G U(iV) may be expanded as 

N 

Det(l - sU) = s k a k {U). (5) 

k=0 



The unitarity of U implies a "self-inversive" property for the secular coefficients ||18|| : 

a N - k {U) = Det(-U)a k (U) . 
Each coefficient a k may be obtained from the traces {U = Till 1 } by 



I fc-i 

i=i 



k-li 



tl 


h 


h ■ 


tk 


1 


h 


t 2 ■ 


tk-i 





2 




tk-2 







k-1 h 

Because this dependence is highly non-linear, the secular coefficients inherit non-Gaussian 
distributions in the RMT ensembles ||. However, to compute the ensemble averages of 
Qu(l) one on ly needs to know their variances, since 

N N/2 

^(l) = J2l k ~ N/2 \a k \ 2 = J2(l k - N/2 + l N/2 ~ k )\a k \ 2 • (6) 

fc=0 fc=0 

For the Poisson and the CUE ensemble of random matrices, these variances were computed 
in 0, and have the following large- iV asymptotics: 

/N \ /I ia\ 

'Poisson \ U I ' \l afc l /( 



;ki 2 > 



'CUE 



Poisson 



CUE 



N 



sin(x/2) 
x/2 



(7) 
(8) 
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In 0, a semiclassical estimation of the \ak\ 2 was given for integrable and chaotic quantum 
maps. The authors used the explicit expression in terms of the traces ti, and estimated 
the latter by the Gutzwiller trace formula ([!]). They made a generalized diagonal ap- 
proximation treating the traces t\ as statistically independent variables. To obtain the 
correlation function, one has to estimate the |afc| 2 (and hence the %) up to times k < iV/2, 
where the Gutzwiller formula is non-rigorous. 



2.2 Representation-theoretic content of Qjj 

We now introduce a more group-theoretic expression for the correlation function. Instead 
of performing the expansion (||), we will express Sluil) as a character in a certain irre- 
ducible representation of XJ(2N), which is best described using the physical language of 
fermions. 

Let JF/v be the Fock space for N types of fermions f { , fj. In mathematics is known 
as the spinor representation space of the group Spin(2iV). Then, for any N x N unitary 
matrix U, 

N 

Det(l -U) = Tr^(-l)^ exp £ //(log U)^ . 

The exponential on the right-hand side can be shown to be well-defined in spite of the 
multi-valuedness of log U. To account for both determinants, we use 2N fermions, whose 
creation operators are denoted by /j. • and f_j , j = 1, . . . , N. The integration over <p in 

the integral (|J) projects on the subspace T = Ker(F + — where F± = J2i f±if±i are 
the number operators for the two types of fermion. The correlation function reads 

N 

n F ( 7 ) = Tr^ 7 ^-^ 2 exp ^ (log UUflj +j - f ,) . (9) 

The operator under the trace belongs to an irreducible representation R of the group 
U(2N), realized on the space which has dimension ( 2 ^). This representation may be 

a b 

defined through its Lie algebra version: any skew-hermitian 2N x 2N matrix X 



is represented by the operator 



d 



N 

R(X) = J2 «, fUf+s + h - f+if-j + •/' >f-J + '*>, f-if-j ■ HO) 

By exponentiating, i?(expX) = expi?(X), we obtain a U(2iV)-representation, which we 
still denote by R. The correlation function ^(7) for any N x N unitary matrix U may 
be recast as a character in this representation: 

^(7) = 7 - Af / 2 Det([/)- 1 Tr R(TV) (11) 
where U d = f g ° ) , r = f ^ J) e U(2iV) . 
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As it stands, the construction assumes 7 = e l6> G U(l). It can also be used for other 
values of 7, since R naturally extends to a representation of GL(2iV, C). In the following, 
matrices in bold print will always be of size 2N x 2N. 

The assignment U 1— ► U embeds U(iV) into U(2iV). By this embedding, R restricts 
to a reducible representation of U(iV) on JF, which we simply denote by R(U). To ex- 
press the correlation function, we may also consider the U(iV)-representation R_i(U) == 
det(f/)- 1 J R(f/). 

In the next section, we decompose R(U) (or equivalently R_i(U)) into irreducible 
representations (irreps) of XJ(N), thus expressing the correlation function Qu{l) as a sum 
of U(iV)-characters. 

2.3 Qjj as a sum of XJ(N)— characters 

The crucial mathematical tool to use is the dual pair structure jHJ. The subalgebra 
{X N <g> la I G u(iV)} of u(2iV) commutes with the subalgebra {Ijv <8> x 2 \ x 2 G u(2)}, 
and each is the commutant of the other inside u(2N): they are said to form a dual pair. 
This means that for all U G U(N), the operator R(U) commutes with the set 

= f+tf-i ' 
% 

i 

J = F + + F^-N. 

The operators J , J 1 * and generate an su(2) algebra. The equation J R(U) = R(U)J 
implies that R(U) conserves the total number of particles and hence acts inside the sub- 
spaces T v = T n Ker(F + + F_ - 2p). 

The dual pair structure provides us with a prescription [TSJ to decompose R{U). Inside 
the reduced Fock space JF, we consider the subspace of lowest SU(2) weights, °JF = 

T fl Ker J|, and expand it according to its particle content: °J-" P = f °JF PI Classical 
results of invariant theory, due mostly to H. Weyl [|34] and succinctly summarized by 
R. Howe [0, amount to the following statements: 

• The operator R(U) acts inside each space °J rp , through a certain irrep p p {U) of 
U(iV). Equivalently, R-i(U) acts on this space through p p (U) = Det(£7) -1 p p {U) . 
Furthermore, two irreps p p and p p < are inequivalent if p 7^ p'. 

• The image of °J rp under (Jf) fc is the space k J r P +k c jF p+fe which is either trivial (if 
k > N — 2p) or carries the irrep p p (if k < N — 2p) . The operators Jj-, Jj_ , Jo ac t on 
this tower of spaces according to the su(2)-irrep of dimension N — 2p + 1. 

• The direct sum of these towers exhausts T . 

We summarize these statements in the following diagram. All entries in a given row 
are subspaces containing the same number of fermions; all entries in a given column (or 
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tower) carry the same U(iV)-irrep. We only show the case where N is an even integer 
(the odd-iV case being very similar): 

jrN _ N-pN 



J7N-1 _ N-ljrN-1 N-2<pN-l 



(12) 



.F 1 = 1 F 1 © 

T J T 

JF° = °JF° 

[/ (AT) - irreps : p pi p 2 ■ ■ ■ Pn/2 

The leftmost tower on the right-hand side carries the trivial u(iV)-irrep, so all spaces 
pjtp = (Jj.)p are one-dimensional. 

Each irrep p p (or p p ) may be described by a Young diagram. p p mixes the action of U 
on p fermions /+ with the action of U on p fermions /_. Owing to antisymmetrization, 
it corresponds to the diagram with p rows of length 2 followed by N — 2p rows of length 
one: 

p p (U N ) = Det(U N )p p (U N ) = Uf lN " P] . 

In view of the above diagram, the dimensions of the representation spaces fc jF fc+p follow 
immediately from those of the spaces T v : 



dim( fc ^ fe+p ) = dim(^ p ) - dim^- 1 ) 

(13) 



N\ 2 ( N x " 



p J \p — 1 

By doing the sum over each su(2)-multiplet we can now express the correlation function 
(|ITD in terms of the irreps p p . 

_ N/2+l-p 

n F ( 7 ) = £ Trp p (U) 1 — I , (14) 

z — ' 1 — 7 

p=0 ' 



or, making the substitution 7 = e lx / N , 

N/2 

Tr Pv (U) — ^ 

sin 



^ S1I1 (2lv) 



For large values of iV we may replace the denominator sin(x/2iV) by x/2N. A quick 
comparison shows that this decomposition is actually equivalent to the pedestrian ex- 
pansion (|(|) written down in Section [2.1| . The squared coefficients |afc| 2 now acquire a 
representation-theoretic meaning: 

p 

Vp < N/2 : \a p (U)\ 2 = Tr^ R^{U) = ^Trp fc (E7) , (15) 

k=0 

or equivalently, 

Tr p p (U) = \a p (U)\ 2 -\a p ^(U)\ 2 . 

As it stands, the decomposition into irreducibles (|14|) is not very informative if one 
takes for U the matrix of a quantum map. We have no way a priori to estimate the 
character Tr p p (Un) from semiclassical information, except by using the relationship, via 
the |afc| 2 , to the original traces Tr([/jy-), as was done in [[/J. This decomposition will, 
however, allow us to obtain rigorous results when adopting a U(iV)-isotropic averaging 
centered around Un (see Section |]). 

2.4 Qjj as a coherent— state integral 

Instead of decomposing the character Tri?_i(£7) into irreducibles, we can rewrite it as an 
integral over the symmetric space M.^ — U(2N)/U(N) x U(iV). This integral can be 
interpreted as a variant of the non-linear sigma model used in [pl| to study the spectral 
statistics of quantum chaotic Hamiltonians on infinite-dimensional Hilbert spaces. In our 
case the integral representation is exact, and is well-defined mathematically. 

To write the character flu(l) as an integral, one uses the coherent states i?(g)|0), 
where |0) is the vacuum of T and g any matrix in U(2iV). These coherent states provide 
a resolution of unity on JF, i.e. they can be combined to build the orthogonal projector 
on T ', as: 

p^ d = f J dg Rigmmig)- 1 , 

U(2AT) 

where the Haar measure dg has to be suitably normalized. Let be the block-diagonal 
subgroup XJ(N) x XJ(N) of U(2JV). Then for all h G H N , the states i2(g)|0) and i?(gh)|0) 
only differ by a phase factor. Therefore, it suffices to integrate over the equivalence classes 
in U(2iV) modulo H N : 

Pt= J d[g] H R(g)\0}(0\R(g)-\ 

V(2N)/H N 

It is convenient to represent the if ^equivalence classes (i.e. the points on Mn) by 
2N x 2N matrices. To each g 6 U(2iV) one associates Q s = gS 3 g _1 , where S 3 = I N ® o" 3 . 
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The set of all these matrices Q is isomorphic to M.n- It is the set of all Hermitian 
matrices with two eigenvalues, +1 and —1, each with multiplicity N. This non-linear set 
of matrices is naturally equipped with U(2iV)-invariant symplectic structure and metric 
(and therefore an invariant measure dQ). 

The matrix elements Qij are not all independent, and for practical calculations we 
need to introduce a bona fide coordinate system on M. N . If we denote by Q12, Q22 
the two N x N blocks in the right half of the matrix Q, the entries of the complex 
matrix Z = Qi2(Q22 — I) -1 are good coordinates on the open subset of M.^ where 
{Q22 — 1) is invertible. Geometrically, these N x N complex coordinates represent a 
certain stereographic mapping of M.n onto C NxN . The matrix Z corresponding to a 
point Q g can be extracted from the Gaussian decomposition of g: 



1 Z\ (A 
1 [C D 



(16) 



These complex coordinates also provide a simple definition of the coherent states. 
Indeed, i?(g)|0) is co-linear with 

l^) d =exp{^/|.^/! J }|0)=exp{ J R^ f )}|0> = J2 (J \ ) |0> . (17) 

As it stands, \Z) is not normalized, but has the following properties: 

• The overlap between two coherent states reads (Z\Z') = Det(l + Z^ Z'). In particular, 
the norm of \Z) is Det(l + Z^Z) 1 ' 2 . 

• The resolution of unity takes the form 

Pt= f d» N (Z,tf) 'f_f ' (18) 
Jcnxn Det(l + Z^Z) 

where the measure dfi N (Z, Z^) = Cn X Det(l + Z^Z)~ 2N nfj=i d^Zy/ir is the ex- 
pression for dQ in the coordinates Z^. The value of the normalization factor Cn is 
given at the end of Appendix B. 

• The group U(2iV) acts on these coherent states as follows: 



R 2) \Z) = Det(CZ + D) \{AZ + B)(CZ + D)~ l ) 
The resolution of unity allows to write the character ([TTD as 

^(7) = 7~ iV/2 Det(t/)- 1 ^ d[ S ] H (OlRigy'R ^ °)i2(g)|0> 



(19) 



'Mn 

C NxN M N 
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This expression is the central result of the current section. It is an exact formula, which 
parallels the "ballistic" non-linear sigma model derived in [|TjJ for Hamiltonian systems 
with an infinite-dimensional Hilbert space. In our finite-dimensional framework, the non- 
local field Q(q', q) of 4 x 4 supermatrices on configuration space is replaced by a "lattice 
field" Qiaj/3 of 2 x 2 matrices (with elements indexed by a,/9) depending on two discrete 
positions i, j. The "effective action" of the present model 

5( 7 , U, Q) = -Tr{log(l + jZ^UZU- 1 ) - log(l + ZZ*)} + y log 7 , (21) 

can be presented []22| in the form 

S(e ix / N ,U,Q) = -Tr]og[coak{H XtU )- ^h{H X)U )Q] , 

with H x u = f — H — log U . 
AN 2 

In [EI3L this action was further transformed, using the Wigner representation of wave 



functions, to obtain the same ballistic non-linear sigma model as in [13|. We will not 



perform these steps, which require some further approximations, but rather try to estimate 
the integral with the above (purely quantum) effective action. 



3 Saddle— point analysis of the action 5(7, U, Q 

To estimate the field integral of their non-linear sigma model, the authors in |Tj| expand 
the effective action around two critical points (usually referred to as saddle points in this 
context). Since there is no large parameter in front of this action, a leading-order saddle- 
point expansion - see Eq. M) - is a not justified mathematically a priori. In the present 



section we explicitly compute this expansion for the action (21) and compare it to the 
results of |L3| and the exact correlation function. 



The saddle points are determined by requiring the variation of the action to be zero. 
In the absence of a large parameter, one first needs to understand exactly which action 
to vary. This point is not entirely obvious: one might be tempted to lift (part of) the 
denominator Det(l + ZZ^)~ 2N of the measure dfi^iZ, Z') into the exponent; this modi- 
fication of the effective action would yield a different saddle-point expansion. However, 
the requirement of coordinate invariance tells us to keep the U(2iV)-invariant measure dQ 
as it is, forbidding such manipulations. With this convention the saddle-point expansion 
of 5(7, £7, Q) will turn out to yield the exact ^-dependence for VLjjij). In particular, the 
problem of "unphysical zero modes" occurring in |13], |22| is resolved. 

We now describe the saddle-point analysis of 5(7, U, Q) in some detail. We first use 
the fact that the action is invariant under simultaneous rotations of both U and Q: 

5( 7 , U, Q) = 5( 7 , VUV~\ VQV- 1 ), (22) 

where we used the shorthand notation V = V ®l2, for V G U(iV). Such a ^-rotation of Q 
is an isometry of the Riemannian manifold M.^ and leaves the measure dQ invariant. It 
therefore suffices to study the simpler situation where U is diagonal. U = D = diag(e ie 7 
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One sees from formula (PH ) that the point Z = (or equivalently, Q = S3) is a saddle 
point, and the quadratic approximation to S for small Z reads: 



S(j,D,Z) « |A^log7-Tr( 7 Z t J DZ J D- 1 -Z t Z) 

TV 

« iiVlog 7 +^|%| 2 (1-7^) • 

m=i 

This saddle point is the only one on M.^ which is located at a finite Z. It is sometimes 
called the "perturbative" saddle point in the physics literature. For a generic matrix U, 
there are N directions Zjj that have a coefficient (1 — 7) ~ —ix/N; these directions are 
called "zero modes" [jT3| , because their coefficient vanishes as x — > 0. Doing the integral 
in this quadratic approximation around Z = yields 

Quil) h = ° N (1-7)^^,(1-7^)) ' (23) 

We chose to separate the zero mode contributions from the others. 

The existence of a second saddle point was pointed out (in the context of the diffusive 
non-linear sigma model) in [ 2Sfl . It may be exhibited through the change of variable 



Z' = l/Z, which amounts to switching to the stereographic projection of M.^ from the 
antipodal point. In terms of the new variable Z', the integrand reads 

Det(l + 1 - 1 Z' ] UZ'U- r ) 



Det(l + Z n Z') 

so it has the same structure as the original integrand, but for an additional prefactor 7^ 
and the replacement 7 — > 7 -1 in the determinant. Quadratic expansion around Z' = 
(or, equivalently, around Q = — S3) yields 

a "M\-* = ° N (l-7- 1 ) JV n^(l-7- 1 e^)) • (24) 
These two saddle points Q = ±£3 (we call them "standard" ) are the only ones taken into 



account in the treatment of the ballistic non-linear sigma model in [|13|, |22j . The problem 
with this approximation is that, in the limit 7 — ► 1, the sum of the two contributions 
Eq. (|2~3D and (pi]) diverges at least as strongly as 1/(1 — -y) N ~ l , whereas the exact corre- 
lation function is bounded. This phenomenon was attributed to the iV — 1 "unphysical" 
zero modes appearing at each saddle point (as opposed to the single "ergodic" zero mode 

Zjj). More generally, these contributions become singular each time U and , yU happen 
to have common eigenvalues. 

We will argue below that this problem with zero modes is actually resolved by taking 
into account further saddle points of the effective action. 
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3.1 Weyl character formula 

To identify all saddle points, we return to the expression (p0|) of the integrand. We still 
study the case where U = D is diagonal, and we write TD = diag(7_D, D). 

Let ( be a complex N x N matrix. The point Q s of A4 n is a saddle point of the 
integrand iff the Taylor expansion of (C|i?(g _1 rDg)|C) around ( = contains no term 
linear in ( and (Note that this statement is independent of the choice of representative 
g for Q s .) Moreover, we do not want the integrand to vanish at £ = 0. If we decompose 

the unitary matrix as g _1 rDg — \ K), these conditions read: 



c d 

b = c = , Det(d) ^ . (25) 



N 



a 



This means that the matrix g _1 TDg (for 7 = e iX / N e U(l)) belongs to the subgroup H 
of U(2iV), which in turn allows g to be written as the product of a permutation matrix 
with some element h G H^. By g CT we mean the unitary matrix (g CT )y = 5i,<x(i)j where a 
is a permutation of {1, . . . , 2N}. To each permutation o there corresponds a single point 
Qa = go-Sago- -1 . Moreover, two permutations a, a' lead to the same point if o = o'r 
where r permutes indices separately inside {1, . . . , N} and {N + 1, . . . , 2iV}; this property 
defines a partition of the symmetric group &2N into ( 2 ^) equivalence classes, each one 
corresponding to a saddle point of the integrand. 

These classes are in one-to-one correspondence with the sets S = er({l, . . . , N}), so 
we can write Q a = Qs- Qs is then the diagonal matrix with entries +1 at the positions 
j G S, and —1 at the positions j G S (the complement of S in {1, ... , 2iV}). We partition 
the set S into Si = Sfl {1, . . . , iVjjmd S 2 = Sfl {N + 1, . . . , 2N}. In the following we will 
also use the set 52 = {j — N\j G S2}, and the sets Si and S 2 which are the complements 
in {1, . . . , iV} of Si resp. The point Qs corresponds to the following (coherent) state 
in T: 

i5) d = f J R( gCT )io)=±n/!,n /-ii°>- ( 26 ) 

The matrix g CT admits a Gaussian decomposition ([H]) iff a is in the trivial class, i.e. S = 
{1, . . . , N}, which explains why only the perturbative saddle point Q = S 3 could be 
exhibited from the Z-coordinates. 

We now compute the leading-order contribution from each saddle point Qs- In the 
vicinity of Qs the integrand in (120) takes the values 



(Cli^rDg^iOACiC), 

where the entries of the matrix ( are "small" (( defines a local coordinate system near Qs)- 
We partition the diagonal matrix g~ 1 rDg cr into two halves: g~ 1 rDg CT = diag(A 1; A 2 ). 
The above integrand then reads 

„ Det(l + ^AiCA,- 1 ) 

Det(.D) x - s 2 



Det(l + CtC) 
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Expanding to quadratic order and integrating over (, £t, we obtain from the saddle point 
Qs a contribution similar to fl2~3|): 



' ies 1 ieSi 

jss 2 jes 2 

x n (i-^e-^-^) -1 n (i-e^-^y 1 n . 

Note that the product contains a factor (1 — , -f)~ N+2r , with r = n 5*2). The most 
singular case r = arises for 5*2 = Si, i.e. saddle points of the type YlieSi f+if-il®)' 

In the general case U = VDV^ 1 , the saddle points are the points Qv,s — VQsV -1 , 
and they lead to the same contributions (cf. the covariance of the action and the measure 
dQ). The two standard saddle points Q = ±S 3 are the only ones unaffected by these 
V^-rotations. 

It is illuminating to present the result of the approximation ( |27D in an alternative 
fashion. For that purpose, we denote the non-zero elements of the diagonal matrix TD 
by e 19 ^ [y = 1, . . . , 2N). The sum of contributions fl2~7|) can then be rewritten in the form 

rr _ rr - pK^m-'M 

^■^? iuifr-^) ' (28) 



Save for the pref actor C^, the expression ( p8|) agrees with the result that follows from 
the VFeyZ character formula |23[] for the trace of i?(TD) over JF. In general, this formula 



expresses the character of an element of U(2iV) (more generally, GL(2iV, C)) in some 
representation R as a sum over all permutations a e 6 2 jv (this being the so-called Weyl 
group of U(2iV)). In our case, the terms from the (2N)\ elements of & 2 n m ay be grouped 
into ( 2 ^) classes, according to the equivalence relation described above. Since Weyl's 
formula is an exact result, the expression remains finite in the limit 7 — > 1, which 
means that the singularities 1/(1 — ^) N ~ 2r of the various terms cancel each other. The 
complete sum over saddle-point contributions thus solves the problem of "unphysical zero 
modes", i.e. the divergence problem of the two standard saddle points. 

The mathematical reason behind the "almost exactness" of the leading-order saddle- 
point expansion is as follows. The action of R(g) on coherent states \ Z) may be interpreted 
as the equivariant action of g on the space of holomorphic sections of a certain complex 
line bundle Lr over M.^ @- This equivariant action can be extended to the (infinite- 
dimensional) space of square-integrable differential forms of degree (0,p) on the bundle. 
On the enlarged space, the character becomes a (super)trace, which can still be written 
as an integral over .Mat. Owing to an iV = 2 supersymmetry, the integrand may be 
continuously deformed without changing the value of the integral. In one limit of the 
deformation, one gets Tri?(g); in the other, the integrand localizes at the fixed points of 
g on AIjv, yielding Gaussian integrals around these points. 

It turns out that these fixed points coincide with our Qv,s, an d their (Gaussian) con- 
tributions are equal to (p7| ), save for the prefactor CV As a result, the leading-order 
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saddle-point approximation (for our non-localized integrand) delivers the correct answer 
(omitting the prefactor). In Section [4.4| , we investigate the higher-order terms of the 
expansion at Q = £3 up to two loops: we find that these terms only renormalize the pref- 
actor Cn, without affecting the U- or 7-dependence. We speculate that the (adequately 
resummed) full series yields the exact answer, including the correct normalization. 

To achieve agreement with the Weyl character formula, it was crucial to regard the 
denominator Det(l + Z^Z)~ 2N as part of the measure (as opposed to lifting it into the 
action). Indeed, in order for the mechanism of equivariant localization to take effect, 
the integration measure must be U(2A r )-invariant - a property not enjoyed by the flat 
measure Ylij d 2 Zij without the factor Det(l + Z^Z)- 2N . 



4 Why do we need averaging? 

While the Weyl character formula for fit/ (7) constitutes an exact result, it is of no use 
- at least not as it stands - towards our goal of proving universality of the correlation 
function. This formula relies on the knowledge of the eigenphases e l8i of U, which are not 
given a priori. It does not exhibit the semiclassical features of the quantum map at all. 
On the contrary, it is a "purely quantum" decomposition of the correlation function, a 
complicated reordering of the Fourier decomposition @. 

As was explained in the introduction, it is not conceivable in general that a universal 
result for flu can be obtained without doing some kind of averaging over the matrix U. 
Given the results of the previous section, one might try to perform the averaging term by 
term in the Weyl decomposition, hoping that most of the terms might average to zero. 
Such a hope is quickly discouraged by a look at the expression : aside from having an 
iV th -order singularity at 7 = 1, whose degree increases each time some 7 e 1 ^ -6 ^ crosses 
unity, this contribution to fit/ (7) is strictly positive for real 7 < 1. We know that the 
singularities are artifacts of the Weyl decomposition, as the correlation function flu (7) 
itself is uniformly bounded w.r.t. U and 7. Unfortunately, because of the positivity of 
( p3| ) the singularities can only be removed by reorganizing the entire sum of contributions, 
not by averaging individual terms. 

For this reason, we will adopt a different strategy: we first perform the (1Vn{U) average 
on the integrand of the coherent-state integral, obtaining a new effective action 

e -Sav(7,Q) ^ ( e -s(rf,u,Q) ^ _ ^29) 
We then estimate the resulting Q-integral by performing a saddle-point approximation 



on the action S av (Q) |2T 



A priori, this approximation is no more justified than the one in the previous section, 
as S'av is preceded by no large parameter either. The absence of a large parameter also 
implies that averaging and making the saddle-point approximation are non-commuting 
operations. Therefore, the saddle-point expansion of will yield qualitatively different 
results from the direct expansion for S(j,U,Q). We explained above that averaging the 
Weyl character formula is hopeless for our aims. The other way around (i.e. performing 
the expansion after averaging the action) will prove more interesting. 
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4.1 Where are the critical points of 5 av ? 



For any averaging measure &Vn, the two points Q = ±S 3 remain saddle points of 
5av(7, Q). In the vicinity of E 3 , the integrand expands as 

( Det De!(l 7 + t ^" 1) ) » ((^AiU- Z) - Ztz) = e-^""-^^ 

where Ad [/ ■ Z = f UZU~ l is the adjoint action of {7 on Z. The approximation is valid 
for Z small. For larger values of Z, one should add higher cumulants to the right-hand 
side. However, for the time being we stick to the purely quadratic approximation, and 
carry out the Gaussian integral to obtain 

<fi C7 ( 7 )>| E3 = C i v7" Ar/2 Det(l-7<Adf/»" 1 . (30) 

When the averaging is absent (that is, (IVn is a Dirac 5-measure at U), we recover the 
contribution ( p3|) . The saddle point Q = — S 3 yields the same result, with 7 — ► 7 _1 . On 
setting 7 = e iX / N , the sum of contributions becomes 

/ p-is/2 \ 

(ttu{e ix ' N ))\ = 2C N M — rv . (31) 

x v ;/ |s 3 u-s 3 N yDet(l-e™/N (AdU)) J K ' 

In the next section, we examine the possible occurrence of further saddle points of 5 av . 



4.1.1 Searching for other saddle points 



In Section [D] we located the saddle points of the function Q 1— > (Z (Q)\R(TXJ)\Z (Q)) , 
using the action of the group U(27V) on the coherent states \Z). This function may be 
interpreted as the Husimi function (or Q-symbol) of the operator R(TU) acting on J 7 , 
and we denote it by i?R(ru)(Q)- By the same procedure we can obtain the saddle points 
of Hr( s )(Q) for any non-degenerate matrix g e U(27V); in that case, the saddle points 
Q cr it are given in general by finite matrices Z crit and Z\ T - lt , which are solutions of the 
saddle-point equations 

() tf* (g) (Z,Zt) = = J-tf R(g) (Z,Zt) (i,j = l,...,AT). (32) 



dZy ' ' dZ v 

It is useful to extend i7i?( g ) to a function of two independent complex matrices Z, Z* (that 
makes 2N 2 complex variables). The saddle-point equations pose 2N 2 constraints on the 
degrees of freedom Z and Z*, which yields isolated solutions (Zi,Z*), provided that the 
constraints are independent of each other. 

The reality of these solutions (i.e. Z* = (Z^) is due to a symmetry of the operator 
i?(g), which is not conserved if we replace -R(g) by any operator TZ on T. For instance, 
if the representation R is extended to matrices G G GL(27V, C), one can show that the 
saddle points of H R ( G )(Z, Z*) are real iff G is a normal matrix (i.e. GG^ = G*G). We 
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are presently unable to determine the conditions for the saddle points to be real for the 
most general TZ. In any case, the saddle points will be real if 1Z is a Hermitian operator. 
The Husimi function is then real, and Morse theory applies to it. By Morse's theorem 



26(1 , the number of saddle points (which we assume to be isolated) is at least the sum 
of all Betti numbers of A4n, which is ( 2 ^) f2~?fl . This is exactly the number of saddle 
points we found for H^x.){Q) when X is a 2N x 2N Hermitian matrix, so this function 
is what is called a perfect Morse function for M. N . X can be joined to g 6 U(2iV) by 
a continuous path inside the set of non-degenerate normal matrices: this explains why 
Hr(t\j)i although a complex function, still has ( 2 ^) real saddle points. 

Unlike reality, the property that the solutions of (^) are isolated points is robust; 
(Zi, Z*) are the common zeros of 2N 2 polynomials in Z and Z*, so they are stable 
w.r.t. perturbations of the coefficients, as long as the equations do not become degen- 
erate. In Section |3j] the saddle points of H^ru)(Q) were called Qv,s- We now switch 
to such complex coordinates ( that a saddle point Qv,s is situated at ( = = , and 
perturb R(FU) in GL(T) to K= R(T\J) + e5K. Then for e small, H n {Q, C 1 ") will have an 
isolated saddle point at (( € , Q), where both ( € and Q are of order e. Even if it is not real, 
this saddle point will contribute to the integral over M.^: starting from real coordinates 
3?Cij) ^Cy; we can locally deform the contour so as to reach the point 

= (Cm,- + QJA (3Cy) crit = (C« - C,J/2i, 

and we can compute the saddle-point expansion of J Hn(JR,(, around it. 

The averaged integrands we want to consider are all of the type Hn(Q), where 

n= ( dv N (v) -JL- r(tv) , 

Jv(N) DetV 

and (1Vn{V) is a normalized measure on U(iV). If this measure is very strongly peaked 
near a matrix Un, the resulting operator will be a perturbation of i?(rUN)/Det£7jv, so 
the above stability arguments apply: the saddle points are then isolated points near the 
unperturbed ones, and they are "almost real" and hence will lie on the integration contour 
after a slight contour deformation. 

For less concentrated measures dV^iy), the structure of the saddle points can change. 
In Section || we exhibit an averaging scheme for which the saddle points are real but not 
isolated: they form submanifolds of M.n] this is also the case for Hr^ if g is degenerate. 
We do not have a good estimate of the typical "width" of the measure (1Vn(V) above which 
saddle points can coalesce, spread over higher- dimensional sets, or cease to contribute to 
the integral (for instance when they depart too far away from reality). 

In general, we are unable to explicitly locate these extra saddle points, even for the 
relatively narrow averages described in Sections [O] and || consequently, we cannot do 



better than stick to the approximation ([JT]) to describe the correlation function. The 
remaining task then is to investigate the spectrum of the operator ( Ad U) , which depends 
on U and on AVn- 
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4.2 Common spectral features of (AdU) 

The spectrum of ( Ad £/) has a few features that are independent of the averaging scheme. 
Before averaging, the eigenvalue unity occurs in AdU with multiplicity N, correspond- 
ing to the iV-dimensional space spanned by the [/-eigenstate projectors (j = 
1, . . . , N), and the remaining N 2 — iV eigenvalues lie on the unit circle. After averaging, 
only the uniform mode In = Ylj iV'jXV^'l is left with eigenvalue at unity, while all other 
eigenvalues have moved inside the unit disk. As a result, the sum of the contributions 
( Pip stays finite in the limit 7 — > 1. Averaging thus removes the "unphysical zero mode" 
problem associated with the two standard saddle points in Section [| 

More precisely, the large-iV behaviour of y^ui^ 31 )\±% 3 for finite x mostly depends 

on the positions of the eigenvalues of ( Ad U} closest to unity. Within the approximation 
(Pip, these eigenvalues are the relevant dynamical data of the correlation function. 



4.3 Semiclassical averaging 



In a semiclassical averaging scheme around a quantized map U^n was proposed 
as a promising candidate to obtain universal spectral statistics, differentiating between 
integrability vs. chaotic behaviour of the classical map 0. One chooses a finite set of 
Hamiltonian functions Hj, corresponding to Hamiltonian vector fields S#. (j = 1, . . . , r), 
on the classical phase space. These Hamiltonians are quantized on each of the quantum 
Hilbert spaces T~Cn, yielding operators {Hj}, which are represented by Hermitian N x N 
matrices w.r.t. an orthonormal basis of Hn- An ensemble average is then introduced by 

1. composing U^n with the operator exp ^— i ^\ tjHj/h\ , where the "times" tj are real 
numbers; 

2. averaging over the parameters tj in a window around the origin of width e using, for 
instance, the Gaussian weight (e 2 7r)~ r / 2 e~^j V e . 

The width e is taken to be ^.-dependent: e ~ h a ~ N~ a for some 1 > a > 0, so that the 

probability measure for the classical maps exp (^jtj^H^ shrinks to a single point, 

0, in the classical limit iV — > 00. The set of Hamiltonians {Hj} is chosen once and for 
all, and is independent of N and the map 0. The only constraint on this set is that the 
second-order differential operator —A = 52 7 -— if- mus t be elliptic pi . 



As explained in the introduction, this averaging procedure is introduced in order to 
suppress the non-generic spectral statistics of quantum chaotic systems with arithmetic 



symmetries. In this respect we must mention the results obtained in ||10|| , where the 
authors show how non-linear perturbations of quantum cat maps exhibit generic spectral 
statistics, as long as one perturbs in both directions of the 2-dimensional phase space; in 
contradistinction, perturbation in a single direction may leave one arithmetric symmetry 
intact, leading to non-generic quantum spectral statistics. This need for "phase-space- 
isotropy" of the perturbations is very similar to our ellipticity requirement: A is elliptic 
only if the vector fields E H . span the whole tangent space at every point of phase space. 
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Some recent articles have dealt with the spectral analysis of the operator 

( Ad t/) semiclas , and obtained interesting results concerning its largest eigenvalues. For a 
classically chaotic map, these were shown to converge (as N — > oo) to the Ruelle-Pollicott 
resonances of the corresponding Frobenius-Perron operator |30|. These resonances are 
inside the unit circle, which means that ^ Adt/jvy . t has a finite gap between unity 
and the rest of the spectrum, for N — > oo. The huge majority of eigenvalues tend to 
accumulate on the origin. 

These properties allow us to estimate the contribution from the two standard saddle 
points for the case of a quantum chaotic map. To lowest order in 1/N, 

(Qn h - e™/ N )) N ~°° NCn sin{x/2) (33) 

^ M7 j; semiclas±E3 ~ Det±(I _/ Adf/ \ ) x/2 {66) 

v \ / semiclas' ' 

where Det^ means that the determinant is computed after restriction to the traceless 
matrices, i.e. to the subspace orthogonal to the uniform mode Ijv- Apart from the non- 
universal prefactor, the x-dependence agrees with the CUE result (§) in the limit of large 
matrices. 

In the case of an integrable map, the eigenvalues of ( Ad t/jvy behave differently: 

some of them populate more and more densely a few curves which connect the origin to 
some point on the unit circle (including unity). For this reason, one cannot separate unity 
in Det (l — e lx / N ( Ad U) . . ) from the rest of the spectrum. All we can say is that the 

V \ / semiclas/ 1 J 

approximation (|3T| ) does not yield the CUE formula in that case (in general it does not 
yield the Poisson answer either). 



4.3.1 Warning 

One might be tempted to present formula (^) as a "physicist's proof" of a weak univer- 
sality conjecture for quantum chaotic maps. The reason why it is not a proof is clear: 

• As was explained in Section [4.1.1| , there certainly exist other saddle points of (|29|). 
The calculation of their contributions is a difficult task, which we have not yet 
performed. It is far from obvious why these saddle points should be less important 
than Q = ±£3 in the semiclassical averaging scheme. 

• As was emphasized before, there is no large parameter in front of the effective action. 
Without such a parameter, the correction terms of the asymptotic expansion around 
each saddle point are not small, and their neglect in the formula fl33|) seems to be 
unjustified. 

The second worry is addressed in the next subsection. 



4.4 Loop expansion 

We are now going to investigate those corrections to the formula ( |5B| ) that result from 
systematically expanding around the saddle point Q = E 3 . The computations will be 
done up to what is called two-loop order in field-theoretic language. 
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As a first step, we approximate the integrand by taking the ensemble average inside 
the determinant: 

(Det (1 + ^AdU ■ Z) ) scmidas « Det (l + 7 ^ f (Adf/) semiclas • Z) . (34) 

Although e, the "width" of the perturbation, decreases like h a , its effect is strong enough 
to completely modify the spectrum of AdU, even in the semi classical limit. This shows 
that the above approximation is not necessarily valid if we just suppose that the matrices 
Z are bounded (uniformly w.r.t. N) in the operator norm on T~Cn- Using the expansion 

N 

Det(l + A) = 1 + TrA + Tr(A J 'A) , 

i=2 

Eq. fl34|) will hold as long as the terms for j > 2 are small compared to TrA. A sufficient 

condition for that is Tr(|A|) <C 1, where \A\ == V A^A. Upon the replacement A = 
ryZ'AdU ■ Z, this condition will be met if 

N 

Tr(ZtZ) = J2 l%| 2 « !. (35) 

i,i=l 

uniformly w.r.t. N. It would be desirable to better control the error in (|34|) for the larger 
set of matrices Z satisfying (iV-uniformly) \\Z\\c(ft \ < const. 

Taking (^3) for granted, we proceed to the computation of higher loops. To simplify 

the notation we abbreviate T = f 7(Ad[/) sem icias- Next we formally introduce a parameter 
M (which will be reset to unity at the end of the calculation) by making in the integrand 
the replacement 

Det(l + Z^TZ) /Det(l + &TZ) \ M 
Det(l + ZTZ) ~* V Det(l + ZW) ) ' 

A contribution to the perturbative saddle-point expansion is said to be of n-loop order if 
it varies as M~ n relative to the leading-order term. On rescaling the integration variables 
to ( = ZyM and £t = Z^yM, the 1/M expansion of the integrand looks as follows: 

^•* C»y + 1IS ) =C " n ^|e-<-.<( 1 + M-V I + M- /2 + ...) 

where fi and /2 are the one-loop and two-loop terms, respectively, and are given by 

h = |Tr(C f C) 2 - |Tr(C f TC) 2 - 2^^, 

f 2 = -iTr^C) 3 + IM^TQ 3 + | [XrCC+C) 2 - Tr^TC) 2 ] 2 

+ 2A^ 2 (TrC t C) 2 + A^Tr(C f C) 2 - NTi^Q (Tr(C f C) 2 - Tr(C f TC) 2 ) . 
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The Gaussian integral at leading order just yields the result ([ft]). Using standard dia- 
grammatic techniques to do the one-loop integral we find the following expression: 




By the relation (1 — T) 1 = 1 + T(1 — T) 1 these terms combine to yield the simple answer 



N 



C » J Il^e- T < t ( 1 -^/ 1 (C,C t ) = ^M-^Det(l-T)- 1 (-iV 3 ) . 

We see that the dependence of the one-loop contribution on T cancels completely, leaving 
only a constant, —N 3 . This cancellation is not accidental but continues to higher loop 
order. By a lengthy but straightforward calculation, the complete perturbative result up 
to two-loop order can be shown to be 

^ v ' Det M (l + ZtZ) 

£NxN 

= C N M- NxN Det(l - T)- 1 (1 - M~ l N 3 + M~ 2 (\N 6 + ^N 4 - ±N 2 ) + 0(M- 3 )) . 

Again, all the T-dependence has disappeared from the loop correction terms. This is true 
for all M including the case of interest, M — 1. 

The cancellation does not come as a total surprise. The above perturbation expansion, 
whose low-order terms we have computed, is formally identical to the same expansion 
before averaging. The latter is obtained from the former by simply substituting jAdU 
for T = 7(Ad?7) sem i C i as . In the case before averaging we know from |^4| that an index- 
theoretic mechanism (sometime called localization) causes the perturbation expansion 
to be deformable (by an underlying N = 2 supersymmetry) to a harmonic oscillator 
problem (or, equivalently, a Gaussian integral) at Z = 0. The process of deformation to 
the Gaussian limit explains why the dependence on jAdU is exhausted by the leading- 
order term. It leads to the Weyl character formula, which implies that the contribution 
to the character from Z = (or Q = S3) is exactly given by 



Det(l + ZtZ) 



Z=0, all orders 



where the normalization constant Cn has now been replaced by unity. The last fact 
provides the raison d'etre for the iV-dependent terms produced by the loop expansion: 
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their role is to cancel, after proper resummation, the prefactor Cjv- This property does 
not depend on the unitarity of 7 Ad U, so it holds as well after replacing it by its average. 
Thus, after summing all orders of the perturbation expansion, we expect that the saddle 
point Z = contributes to the correlation function as 

< Si "W>| E3 ,.„„ dc „ = 7- A ' /2 Det(I-o^Adt/))-V 

This perturbative result should be used with some care. Although the function 
f(Z, Z^; T) = Det(l + Z^TZ)/ Det(l + Z^Z) is locally well-defined, it does not extend to 
a global smooth function on the manifold M.^ (in particular, this function is NOT the 
Husimi function of an operator on T\ Indeed, setting Z = zG with any invertible matrix 
G and sending z — > 00 always leads to the same point Q = — S 3 on A4n, regardless of 
which matrix G we choose, whereas the limit of f(zG, zG<; T) as z —>■ 00 does depend on 
the choice of G. Thus, the function f(Z, Z'\ T) is not smooth at Q — —S3. 

This singularity reflects the fact that the cumulants neglected by our basic approxi- 
mation ( P^D are small (compared to the terms kept) only for small matrices Z (cf. the 
discussion following Eq. (0)). If Z, Z^ (or some matrix elements thereof) are allowed to 
go to infinity, the approximation clearly loses its validity. To control the error incurred 
near the saddle point Q = — S3, one needs to switch to another scheme, by first changing 
coordinates Z — > 1/Z and Z^ — > 1/Z* and only afterwards repeating the above steps. 
The contribution from Q = — S 3 can then be calculated in the same way as the one for 
Q — S3. The treatment of further saddle points remains an open problem. 

What makes this procedure unsatisfactory is that we are simultaneously working with 
several approximation schemes, each of which is only locally controlled. To localize the 
integral at the saddle points in a mathematically rigorous manner, we would need an 
approximation that is globally well-defined and well-controlled. It is not clear whether 
such an approximation exists, given the stringent requirement that the integrand should 
also have the index-theoretic features that allow localization techniques to be used. 



5 Averaging U over eigenbases 

By its definition (0) as a correlation function of spectral determinants, fit/ (7) is in- 
variant under any change of basis U 1— > V77V -1 , with V an arbitary unitary matrix. 
In the Q-matrix formulation, this invariance is reflected by the relation 5(7, U, Q) = 
^(7, V77V -1 , VQV -1 ). Since the transformation Q 1— > VQV^ 1 has unit Jacobian, we 
may absorb V into the integration variable Q and compute flu (7) by first averaging 
e -s(i,u,Q) over gji ro t a tions U i-> VUV" 1 : 



e -s Vav ( 7 ,t/,Q) 4| f 



Vol U(JV) 



I dV expi-S^VUV- 1 ,®)}, (36) 

JV(N) 



and then integrating e~ Sv ^ v over Q. We saw in Section |3.1| that if the matrix D = 
V^UV is diagonal, then the saddle points of 5(7, U, Q) are situated on the points Qv,s — 
VQ5V -1 . Because the locations of these points explicitly depend on V, we expect that a 
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smoothing mechanism takes place and the divergences of the individual terms in the Weyl 
character formula disappear on averaging over V. In fact, as we will see, the expansion 
obtained by saddle-point analysis of the effective action SVav(7, U,Q) is qualitatively 
quite different from Weyl's formula. 



5.1 Analysis around ±£3 

We first describe Svscv(j,U,Q) near the two saddle points Q = ±S 3 (cf. Section |4.1|) . 
The ^/-averaged adjoint operator ( Ad ) y has a rather simple spectrum: unity is a simple 
eigenvalue (associated with In), and on the remaining (N 2 — 1) -dimensional space the 
operator is proportional to the identity: 

(MU) V = I\ + (1-I\) \™} 2 -~i ■ < 37 ) 

(Pi is the orthogonal projector on In-) We see that (Ad£7) y has a large gap between 
unity and the second eigenvalue, and this gap has the maximal degeneracy. Assuming 
that this degenerate eigenvalue is small (|Tr[/| <C N), we get the following leading-order 
contribution: 

/o 2NC N sin{q;(l/2-a)} , , 

^ u{G )>v|*u-* ~ (1 - a/Nr x ' (38) 

with a = (\TrU\ 2 - 1)/N . (39) 

Within this approximation, the correlation function depends on U = U^^n only through 
the simple quantity |Tr[/| 2 , which can be estimated semiclassically by the Gutzwiller- 
Tabor trace formula (0): typically, a is of order 0(1/ N) for a chaotic map, and of order 
0(1) for an integrable one. 

Notice that, due to the high degeneracy of the second eigenvalue, we do not get in 
general the CUE result (§), although this eigenvalue is far inside the unit circle. This 
shows that, to obtain the CUE result (|33|), we not only need a finite gap in the spectrum 
of (Adt/(, but also a fast accumulation of the eigenvalues to the origin. The precise 

condition on the eigenvalues is J2f =2 j- 2 ^ ^ N. In the present averaging scheme, this 
means a«l. 



5.2 Critical submanifolds 

We need to investigate the possible influence of other saddle points of 5V av (7, U, Q); for 
the present averaging scheme, we will explicitly describe a critical set, which we believe to 
be exhaustive. The effective action possesses the symmetry SVav(Q) = SVav(W<5W -1 ) 
for all W G U(iV). Therefore, the saddle points are grouped into stationary submanifolds, 
each of them invariant under XJ(N). 
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5.2.1 Description of the manifolds 

In Appendix A we prove the following statement: for any initial matrix U and any 7, the 
action (|36|) is stationary at the points Q = WQo-W -1 , for all rotations W G U(iV) and 
any permutation a (see Section |3J] for the notations a, S = Si U S2 etc). Since U(iV) is 

connected, the points Ais = {WQs-W^ 1 | W G U(iV)} form a connected submanifold 
of M N . 

Let t be any permutation among A indices. The set S' = (r(Si), t(S , 2)) is in general 
different from S, and we have Qs 7^ Qs'; however, Qs> G Ais, or equivalently A4s = A^s/. 
Putting p = §S 2 and r = j^S^flS^), we find that Ais contains ( ) (^) ( ~ p ) different points 
Qg/. The manifolds A^s are in one-to-one correspondence with the integers (p,r), and 
their total number is (A/2 + l) 2 for A even, and (N + 1)(N + 3)/4 for A odd (including 
the isolated points ±E 3 in the count). 

For generic U and 7 ^ 1 (genericity means here that the matrix diag(7?7, U) is not 
degenerate), we conjecture that the submanifolds M.( p , r ) exhaust all the critical points of 
the action SVav(7, U, Q). 



5.2.2 Contributions of the manifolds 

The leading-order contribution of each submanifold Ais to the Q-integral is calculated by 
separating the tangent space at into two parts, one parallel and another one transverse 
to M-s- The integrand in the vicinity of Qs then reads (to quadratic order) 

e -Sv,AQs) e -Hess T (X T )+C>(|X T | 3 ) 

where Hessr is the Hessian of SVav around A^s, viewed as a non-degenerate quadratic 
form on the transverse part of tangent space (coordinatized by Xt)- The exact integral 
over A^s and the Gaussian integral over the transverse directions yield the contribution 

^wV^ c ^" N/2+P v^y e_W7 '" (40) 

In Appendix B, we explicitly compute the volumes of the submanifolds Ais = Ai^^y. 

(r(i) • • -r(r)) 2 r(i) ■■■r( P -r) r(i) • • • r(N -p-r) 



VoLM 



(p,r) 



r(i)---r(A) 



For all submanifolds Ais 7^ {±S 3 } (i.e. < p < A), these volumes are A-exponentially 
small. The quantities Hess^ and SVav(QsO depend on U and 7; we are unable to compute 
them in general. What we know for sure is that | e -5vav | < 1, since e~ s ^' u '^ has this 
property. 

For a non-degenerate U and 7 = e lx l N , the Hessian around Ais will possess a single 
eigenvalue that vanishes with x, while all other eigenvalues stay at least of order 0(1). 
This means that the contribution from At 5 goes like 1/x as x — > 0. However, the "particle- 
hole duality" between the submanifolds A4(p, r ) and Ai(N- P ,r) cancels this divergence in 
the sum of their two contributions (as it does for ±E 3 ). 
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As a result, we conjecture that each contribution \^f7(7))|^ uM _ * s ^-uniformly, N- 

exponentially small compared to that from ±E3 for large N, owing to the small volumes 
of M.{p, r )- Since the number of critical sub manifolds grows like N 2 , we deduce that the 
leading-order saddle-point expansion for the action Sy av ({7, 7, Q) can be truncated to 

@D for lar S e N - 



5.3 Averaging over Random Matrix ensembles 

We may go further and average e~ s not only over the conjugates of a fixed matrix U, 
but also over the spectrum {e l6,J }. For instance we can average U over all matrices in 
XJ(N), with a weight corresponding to one of the standard random matrix ensembles 
(Poisson, CUE). The averaged action will be U(iV)-rotation invariant, and its saddle 
points will still lie on the sub manifolds Ais- As a result, the leading-order saddle-point 
(l.o.s.p.) expansion for such ensemble-averaged actions can again be truncated to the 
contribution (|38|) , upon replacing the coefficient a by its average ( a ) ensemble over the 



ensemble considered. 



5.4 Conclusion: no l.o.s.p. expansion for the V— averaged actions 

The contribution fl55| ) depends in a very simple manner on the matrix U, namely only on 
its first trace. This is in contradiction with the fact that a priori, all traces up to Tr(U N ' 2 ) 
enter into Qu(j) (cf. Eq. (0)). By selecting some particular cases, it becomes obvious 
that the l.o.s.p. expansion fl38|) deviates strongly from the exact correlation function. The 
most immediate counterexample is the Poisson ensemble, whose correlation function is 
given in Eq. (|8]). For this ensemble, ( a ) Poisson = lj which yields the CUE result (!) when 
inserted into the formula (|38|) . We are hence forced to abandon the l.o.s.p. expansion for 
the ^/-averaged actions. 

Nevertheless, we hope that this expansion is still meaningful when the averaging over 
U is local in U(N), which is the case for the semiclassical average in Section [O] (but 
not for the ^-average). Hopefully, a local average will still conserve some memory of the 
"localization" property, which entailed the "almost exactness" of the l.o.s.p. expansion 

fa s(>y,u,Q). 

In the next section, we will consider a local averaging scheme different from the semi- 
classical one. It possesses group-theoretic properties, which will allow us to analyse it 
from the character decomposition ( |14|) instead of the coherent-state integral. 



6 Isotropic averaging 

Starting from a fixed matrix U, one may define an isotropic averaging around U, by 
composing U with the N x N unitary matrices e~ lH , weighted by exp(— TrH 2 /4e)dH 
with small e (so that the weight is concentrated at the identity). Isotropy here means 
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that the measure dH is U(iV)-invariant. Note that this in sharp contrast with the semi- 
classical averaging of Section [4.3| , where H was a linear combination of / matrices Hj, 
with / independent of N. In the semiclassical case, the perturbation spanned only a 
/-dimensional submanifold, whereas in the present case the perturbation completely fills 
the iV 2 -dimensional e-ball centered at H = 0. 

One can replace the Gaussian weight by any positive normalized U(iV)-invariant func- 
tion of H. For our purposes, it is convenient to use the heat kernel on XJ(N), i.e. the kernel 
of the regularizing operator exp(— eA), where A is the (positive) Laplace-Beltrami oper- 
ator on U(N). The heat kernel centered on U is defined as follows: 

Ve>0 : -A v K e (V,U) = %-K e (V,U) 

oe 

lim K e {V,U) = Su{V) . 

Owing to the compactness of U(N), the density K e (-, U) for any matrix U converges to the 
uniform density on U(iV) as e — > oo. Switching e from to oo therefore realizes a crossover 
from the Dirac delta measure 5u(-) to the Haar (or CUE) measure. For small values of e, 
the kernel K e (V, U) = k e (VU~ l ) is concentrated around e~ lH = VU^ 1 ~ 1 and is approx- 
imately given by the Gaussian weight introduced above: k e (e~ lH ) ~ exp(— TrH 2 /4e). 

Schur's lemma ensures that A is proportional to the identity on each U(iV)-irreducible 
subspace of L 2 (U(iV)). As a consequence, its action on each representation matrix p p (U) 
of Eq. (|14]) is simply multiplication by a positive factor, called the quadratic Casimir 
invariant, which we denote by p p (A). In formulas, 

dV Pp {V) K e (V, U) d ^ e- A - p p (U) = e"^( A ) p p (U) . 

U(AT) 

The factor p p (A) may be computed from the Young diagram of Pp ; a more direct way is 
to express A in terms of fermionic operators acting on the Fock space T: 

N 

= (iV + l)(F + + F_)-(F 2 +F 2 )-2J T J i . 
Applying this to any element of the subspace °J rp (which carries Pp ) we find 

p p (A) = 2p(N + l-p) . (41) 

On employing the decomposition (0), the heat-kernel averaged correlation function for 

7 = e iX l N takes the form 

^ air, \Z<^ - 2 -2zl 



(n u ( 1 )) e d ^ f e-^ Stub) = J2e- 2ep{N+1 - p) Tr Pp (U) ^MiZ 



N 



p=0 Sm V2N> 
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The effect of the averaging procedure is to damp the large-p traces, which are difficult to 
estimate from the Gutzwiller trace formula. In the above equation the e — > oo behaviour 
is obvious: all traces except the trivial one Trp (?7) = 1 are killed by the exponential, no 
matter what the matrix U is. It is actually not necessary to set e to oo to get the CUE 
correlation. Since the irreps p p are unitary, their traces are bounded by 

\T*Pp{U)\ < dimpp = Trpp(I). 

The dimensions of the p p s are given in Eq. fll3|) ; for finite p, they are bounded by 
dimpp < N 2p . In the limit N, p — ► oo with y = p/N fixed, Stirling's formula yields 

dim( Pp=Ny ) ~ fr^-fr^*"™ , 

2/(1 - y) 

where the function f(y) = —ylogy — (1 — y) log(l — y) increases monotonically from 
/(0) = 0to/(l/2)=log2. 

For any sequence {Un}n£N, if we tune e (possibly varying with N) such that 

e = iVe»l, 

all the terms making a significant contribution to ([42]) satisfy p <C N. The x-dependence 
of all these terms is the same (being given the CUE correlation x~ l sin(a;/2)), so the 
averaged correlation will also have this dependence. Only the prefactor will depend on 
the matrices Un explicitly. If e is increased further to e ^> logiV, the prefactor itself 
becomes universal. 

These statements hold even in the most general case, when the sequence {Un} is com- 
pletely arbitrary. Therefore, to be able to differentiate between integrable and chaotic 
quantum maps, one must tune the "disorder strength" e to smaller values, so that contri- 
butions from the "high" traces Ttp p (Un) start entering into the answer. To recover the 
Poisson behaviour for integrable maps, one actually needs contributions to flUf) coming 
from the whole region p < N/2. 

This puts us in a no-win situation. On the one hand, we should tune e to small enough 
values so that the high traces p ~ Ny (y > 0) survive and Poisson behaviour stands a 
chance to emerge. On the other hand, for a chaotic map we have no control over these 
high traces (we don't for an integrable map either). 

For our purposes, the present averaging scheme is probably "too algebraic" , as opposed 
to the semiclassical average presented in Section |4.3| . To motivate this statement in the 
spirit of Section |4.3| , let us compare the spectra of the operators ( Ad for the two 
schemes: 

• The spectrum of ( AdU) . , qualitatively depends on the nature of the classical 

r \ / semiclas a a i 

dynamics (see Section |4.3| ). It has a finite gap for a chaotic map, whereas eigenvalues 
accumulate near the unit circle for an integrable one. 

• In the isotropic scheme, (AdC/) is decomposed into the irreps Pq(U) © Pi(U). 
Therefore, apart from the single eigenvalue unity, (AdC/) £ has the eigenvalues 
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{e~ 2e e 1 ^* - ^'}, where {e l9 i} are the eigenvalues of U; the eigenvalue e~ 2e is (N — 1)- 
fold degenerate. This spectrum is qualitatively the same for chaotic versus integrable 
systems. 



(Trp p (t/)) Poisson 



6.1 Crossover Poisson-CUE 

We now present an application of the above scheme in the area of random matrices. 
More precisely, we use the isotropic averaging to build a crossover between the Poisson 
and CUE ensembles, and we derive the transitional determinant correlation function that 
interpolates between the formulas (H). This crossover, as well as the method used to 
compute (Of/V can be compared to the GOE — > GUE crossover studied in 0. 

Our crossover is defined as follows. We start from the Poisson ensemble, then convolute 
it with the isotropic (heat kernel) measure of width e: 

( n u(l)) Poissone = [ dP Poisson (U) [ dVK e (V, u)n v ( 7 ) . 

JV(N) JXJ(N) 

For e = 0, this is the Poisson ensemble. In the large-e limit, the second integral converges 
L/-uniformly to the CUE correlation function, so the output (^u('y))p oisson e does too. 

To calculate the correlation function along the crossover, we will use the decomposition 
( f42"D as in the previous section: averaging being a linear operation, we only need to replace 
the characters Tr p p (U) by their Poisson averages (see Eqs. (H) and fllED ): 

'N\ ( N 
.Pj \P ~ 1, 

The asymptotics of these traces in the regime p, N — > oo with y = p/N fixed, again follows 
easily from Stirling's formula: 

<-&/v(tf)>M„ ~ ^ N r 1/2 ^=i e " m ■ (43) 

The sum over characters therefore approaches the following integral (as iV — ► oo): 

(n p (e^)) p . ^ [ 1/2 dy f{V) ^ X -V X ) e N(f(yy2ey(l-y)) {u) 

\ UK ;/ Poisson , e v ^nJ q y y/y(l-y) x V ; 

In the limit iV — > oo, this integral is determined by the saddle points (rather, the maxi- 
mum) of f e (y) = f f(y) — 2ey(l — y) on [0, 1/2]. Three cases have to be distinguished: 

• If e < 1, the boundary point y — 1/2 is a maximum of f e and is the only critical point 
on [0, 1/2]. Because of the vanishing of the integrand at y — 1/2, the saddle-point 
analysis requires some care. On scales of order O(N ) the result turns out to be 
independent of x: 

which shows that the Poisson result 2 N is retrieved in the limit e —>■ 0. The correlation 
functions starts depending on x on scales of order x ~ 0{N 1 / 2 ). 
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• If e > 1, the maximum of f £ is situated at the point y £ G (0, 1/2) which solves the 
transcendental equation f' e {y) = 0. The correlation function depends oni~ O(N ) 
as 

W"))^ « . ,45) 

The flat correlation function has been replaced by an oscillatory function, with the 
period of oscillation being controlled by the "frequency shift" y £ . When e becomes 
large, the shift vanishes as y e ~ e~ 2e , so the CUE correlation function is retrieved. 

• If e — 1, the correlation function is "critical" (in the sense of a phase transition), as 
the two points y = 1/2 and y e coalesce for e — > 1 to form a degenerate critical point. 
In this case the correlation function varies on scales x ~ 0{N 1 ^). 



7 Conclusions 



In this paper we have adapted the NLcrM approach introduced in ]T3|, [14] to the framework 
of quantized maps on a Hilbert space of dimension N ~ h . We focused on the spectral 
determinant correlation function fluij) instead of the pair correlation function, thereby 
obviating the need to introduce supersymmetry; we obtained an exact expression for the 
correlation function as an ordinary integral over a iV 2 -dimensional complex manifold. 
Because the manifold is compact and the integrand uniformly bounded, no regularization 
needs to be introduced (unlike in [[13]). 

To estimate this integral we expand the integrand around its saddle points, first re- 
stricting ourselves to the leading-order perturbative expansion around each point. Owing 
to the absence of a large parameter in front of the effective action, this approximation is 
uncontrolled, and the connection between its output and the exact value of the integral 
seems fortuitous at best. 

Yet, for any matrix U G XJ(N), we find that the result from lowest-order saddle-point 
expansion of the effective action ^(7, U, Q) coincides with the exact correlation function, 
up to a global prefactor: 

^U(J)U = C N fi C /( 7 )exact • (46) 

l.o.s.p. exp. 

This remarkable coincidence is linked to a cancellation property of the higher-order terms 
of the perturbation expansion, which modify only the prefactor, and is explained by the 
group-theoretic structure of the integrand and the Weyl character formula. Unfortunately, 
the expansion is of no use for estimating the correlation function of quantized maps in 
the semi classical limit. 

We argue that a decent semiclassical estimate of the correlation function Qu(j) can 
only be reached if one takes an average over a set of unitary matrices in the vicinity of U. 
To estimate this averaged correlation, we first average the integrand e~ s ^' u '^ over U, and 
then perform the saddle-point expansion of the output. Because averaging and saddle- 
point expansion are operations that do not commute, this procedure yields an expansion 
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different from that of the "individual" action. At the same time, averaging a priori breaks 
the group-theoretic structure, and with it the exactness (modulo prefactor) of the leading- 
order saddle-point expansion. Moreover, the explicit computation of saddle points and 
their contributions is in general a non-trivial task for a general averaging scheme. 

We have been able to locate the complete set of critical points only for a certain type 
of average, namely averaging over all bases of Hilbert space. This produces a U(iV)- 
invariant effective action, the critical points of which are grouped into submanifolds, 
and are independent of the matrix U we started from (as long as its spectrum is non- 
degenerate). Two of these submanifolds are isolated points; we conjectured that the 
contributions from these two "standard" saddle points, which can be computed explicitly, 
always dominate the leading-order saddle-point expansion. 

The contributions from these two points are unfortunately "too simple" to constitute 
a good approximation of the correlation function, except in some exceptional cases, which 
we do not truly understand. If we average over U G U(iV) with the Poisson measure, the 
saddle-point result strongly differs from the exact one. We are thus led to conclude that 
the leading-order saddle-point expansion of rotation-averaged effective actions does not 
yield a good estimate of the full integral. 

What happens in the case of a local average, i.e. when the weight of the probability 
measure is concentrated near the quantized map Un, is unclear. For one thing, we are 
only able to exhibit the two standard saddle points of the averaged action, but there 
surely exist many more. 

In the case of the "semiclassical" averaging scheme, expansion around these saddle 
points yields results similar to those obtained in fl3|| , except that the "resonances" we 



identify are eigenvalues of a quantum operator. Yet, these resonances for large N seem 
related to the (classical) Ruelle-Pollicott resonances [^, [5SJ, in particular they indicate 
whether the classical dynamics is chaotic or integrable. 

To connect these resonances with the determinant correlation function on a rigorous 
footing, we need two non-trivial assumptions to be fulfilled. First, we must assume that 
the leading-order saddle-point expansion of the (local average) S , sem icias(7, Un, Q) makes 
sense, i.e. gives a good approximation of the exact result; the two-loop calculation around 
±£ 3 in Section £0| seems to support this assumption. Second, hindered by our inability 
to compute the contributions from further saddle points, we are forced to assume that 
the full expansion can be truncated to the two standard saddle points, or at least that 
this truncation provides a reasonable approximation. We presently see no way to prove 
these assumptions. 



Appendices 

A. Proof of criticality of the submanifolds A4s 

To prove that the l^-averaged integrand e~ s v^(-y,u,Q) j g stationary on the submanifolds 
A4s C M-n, we employ the coherent-state formulation of the Q-integral. The point Qs 
corresponds to the state \S) = -R(g CT )|0), and the points in a neighbourhood of Qs may 
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be parametrized as R(g a )\(), where ( runs through the (small) N x N matrices and |£) 
is the corresponding coherent state. The permutation a G is chosen in such a way 
as to interchange the sets S\ and S% = S 2 + N, and to keep Si and S 2 = S 2 + N fixed. 

We write the 2N x 2N matrix g o r 1 rUg cr in the block form , and first compute 

the value of the integrand in the vicinity of Qs before averaging: 

(C|^(g- 1 rug CT )|c> 



(CIO 



Bet(D) (1 + Trp- 1 ^ + BD^) + 0(\(\ 2 )) . 



Then we perform the U-average on U (recall that V G U(iV) acts on U G U(iV) by con- 
jugation: U 1— > VUV~ r ), and study its output on the right-hand side of the above equa- 
tion. To first order in ( and £t, we need the averages (Det(.D)) y , (D~ l C Det(_D)) y and 
{BD- l T>et(D)) y . By decomposing the sets {1, . . . , N} = S1US1 and {iV + 1,... , 2iV} = 
S2 U S 2 , the N x N matrices B,C,D may be written in block form: 

lU Sl§1 \ n _(lU- SlSl \ n _/ 7 % 



where each entry U ss i is a matrix of size Jjs x JJs', whose indices take values in the sets 
s, s'. Thus the ^/-averaged coefficients of the term linear in ( are the following matrix 
elements: 



(Det^jDet^sJ (u^gU^) ) v , (Det^jDet^) (u^Us 2 s 2 ) lm ) 



v ' 



where we have displayed only the dependence on U (and omitted the 7-dependence) . We 
now use the invariance of the Haar measure dV under (left) multiplication by any unitary 
matrix and any diagonal unitary matrix 5 = diag(5i, . . . , 5n) in particular. Under such 
a left translation, the above matrix elements acquire extra factors 5i/5k (resp. 5i/S m ). 
Hence 

<Det( J D)( J D- 1 C) ife ) y = (DetpXZr 1 ^)^/^ for any 8 h 8 k . 

Since i G Si and k G Si (resp. I G S 2 and m G S 2 ) are never equal and the ratio 5i/Sk 
may take any value in U(l), we conclude 

(Det(D)(D- 1 C) tk ) v = 0. 

By the same reasoning, the terms linear in Q vanish after ^/-averaging. 

We have thus shown that the point Qs on M.^ is a critical point of the V^-averaged 
action SVav, Eq. fl36|). By the U(iV)-invariance of SVav, it follows that the whole subman- 
ifold M.s is critical for the U-averaged action, no matter what U is. 
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B. Volumes of the critical submanifolds 



We treat the general case with (jSi = iS-2 = N—p, (jSi = ({£2 = £>, tK'SiFlS^) — (K^in*^) — 
r, and to build Qs we use the same permutation a as in the previous appendix. 

The manifold M. s is given by the set of states {i?(Vg CT )|0) | V E XJ(N)}. These states 
may be written (up to normalization) in the form R{g a )\Cv) where the coherent state \(y) 
is determined by the matrix 



according to Eqs. (|T6| ) and ([17]). The block structure of this matrix derives from the sets 
(Si, S 2 ) vertically, and (5*1, 62) horizontally. 

When V runs through XJ(N), the upper- left matrix Q 1 ^ takes all possible values in 
(^(N-p)xp^ rpj^g ma ^ r j x £(2) j g no ^. independent of . For a fixed Q 1 ', we need to identify 
the remaining degrees of freedom in which is quite easy to do if Q 1 ' = 0, i.e. if V 
has the structure V = diag(Vs 1 5 1 , Vg^). The matrices V S2 § 2 and V§ 2 § 2 in this case block 
decompose as 

V, , = ° ) V,s= ° 



where the index 12 refers to the set S\ D S%, etc. The degrees of freedom of the lower- 
right part of (v are thus two matrices, C^ 11 ) = ^212^1312 e C rx( - N ~ p ~ r \ and = f 
^12 12 ^li 12 e C^~ r ^ xr . They are independent of each other, and take all possible values 
in their respective vector spaces. Since the subgroup U(iV — p) x U(p) of U(iV) acts 
transitively on the submanifold Q 1 ' = of /As, there exists a natural choice of invariant 
measure on that submanifold. It has the factorized form 



: (n)t C (li))iV-p-Q d 2 C (ii) /7r x Det(l + C (II)t C (n) ) P n rf2 Cf ) / 



The matrix Q 1 ' parametrizes a coset space U(iV)/U(iV — p) x U(p), with the corresponding 
invariant measure being Det(l + £W ('( 1 )) Ar FJ . d 2 Qj /it. By group invariance arguments, 
the volume element of .Ms (normalized so that it agrees with the Riemannian measure 
inherited from the Riemannian manifold /An) is the product of the measures for Q 1 ', 
Q 11 , and above. Using this fact and the result |H| 

def f tttt^^i , ^t^-»-» r(i)...r(n)r(i)..-r(m) 



J( m , n ) d ^ f / TT TT ^Det(l + Z^Z)- n 

J C mXn f_ A I* 7T 



we obtain the volume of .Ms: 

VoLM 5 = I(jp, N—p) I(r, p — r) I(r, N — p — r) 
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A similar integral yields the normalization factor of the measure d\i^{Z, Z^) on the 
full manifold Mn- 

i_ _ f fr t w _ x _ r(2)-r(;v + i) 
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